Non-Fermi Liquid behavior at the Orbital Ordering Quantum Critical Point in the 

Two-Orbital Model 
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The critical behavior of a two-orbital model with degenerate d xz and d yz orbitals is investigated by 
multidimensional bosonization. We find that the corresponding bosonic theory has an overdamped 
collective mode with dynamical exponent z = 3, which appears to be a general feature of a two- 
orbital model and becomes the dominant fluctuation in the vicinity of the orbital-ordering quantum 
critical point. Since the very existence of this z — 3 overdamped collective mode induces non-Fermi 
liquid behavior near the quantum critical point, we conclude that a two-orbital model generally 
has a sizable area in the phase diagram showing non-Fermi liquid behavior. Furthermore, we show 
that the bosonic theory resembles the continuous model near the d-wave Pomeranchuk instability, 
suggesting that orbital order in a two-orbital model is identical to nematic order in a continuous 
model. Our results can be applied to systems with degenerate d. 
superconductors and bilayer strontium ruthenates Sr3Ru207. 



and d yz orbitals such as iron-based 



Introduction - A key puzzle with the iron-pnictide su- 
perconductors is one of size: the 0.3% change in the lat- 
tice constant at the structural transition [1] is not com- 
mensurate with the subsequent massive reorganization 
in the electronic system as evidenced most notably by a 
transport anisotropy [2] that can exceed a factor of two. 
Similar incommensurate changes are also seen in the Hall 
and Seebeck coefficients [1] as well as an enhanced tunnel- 
ing signal at zero-bias in point-contact spectroscopy [3]. 
While on theoretical grounds such physics might be ac- 
countable for in the spin sector alone, the pnictides con- 
tain an additional orbital degree of freedom which, when 
present, has been used successfully to explain the dis- 
crepancy between the electron transport and the tiny lat- 
tice distortion in systems such as the manganites and the 
ruthenates [4-7]. The reason is simple. Orbital degrees of 
freedom are part of the spatial symmetry, not an internal 
symmetry possessed by the spin sector. Relying on the 
spin to generate transport anomalies would rest then on 
the magnitude of the spin-orbit effect on the atom that 
possesses the orbital degree of freedom, Fe in the case of 
the pnictides. The spin-orbit coupling in Fe, however, is 
not sufficient to give rise to such transport anisotropics. 
The same is true in the ruthenates and the manganites. 
Further, as is well known from the manganites, coupling 
fluctuating spins with the lattice can only yield modest 
changes in the transport properties [8]. 

We now know from the crucial work of Kugcl and 
Khomskii [9] in the context of multi-orbital Mott sys- 
tems, that orbital degrees can acquire dynamics and 
hence can order in a manner identical to SU(2) spins. 
Orbital ordering, or equivalently orbital polarization, al- 
though driven by a small lattice distortion, can yield 
sizable transport effects in the electronic sector. Based 
on the success of the orbital ordering program in multi- 
orbital systems such as the manganites and the ruthen- 
ates, we[10, 11] as well as others[12-15] have advocated 
that similar physics applies to the pnictdes, though not 
Mott insulators exhibit many of the characteristics of bad 
metals. In the pnictides, as a result of the C4 symmetry 



in the high-temperature phase, the d yz , and d xy orbitals 
are degenerate. Unequal occupancy of the two obritals 
lowers the lattice symmetry to C2 and sizable rearrange- 
ments obtain in the electronic sector consistent with ex- 
periment. For example, we have shown[16] using the 
random-phase approximation that orbital fluctuations 
between the d xz and d yz orbitals in a five-band model[17] 
for the pnictides can lead to a break-down of perturba- 
tion theory and drive an instability to a non- Fermi liquid 
state. Since the hunt for non-Fermi liquid states is in its 
infancy, it would be advantageous to establish the onset 
of non-Fermi liquid behaviour from the simplest model 
possible and from a non-perturbative method. Such a 
demonstration would help establish what deviations from 
a one-band model lead to an effective breakdown of the 
Landau quasi-particle picture. It should be stressed that 
our result is different from [18] because the dispersions for 
d x 

are present 



and dy Z orbitals do not intersect when interactions 



In this paper, we approach the problem of the emer- 
gence of non- Fermi liquid states of matter using multidi- 
mensional bosonization[19-22]. Since we are after univer- 
sal physics, rather than the starting from the complexity 
of a five-band model, we focus just on a two-band model 
with degenerate d xz and d yz orbitals to see if orbital fluc- 
tuations can give rise to non-Fermi liquid behaviour. We 
demonstrate here such a model can be solved analytically 
and produces the desired result of a non-Fermi liquid 
state of matter. We establish that at the C4 symmetry 
breaking quantum critical point, a z = 3 overdamped 
collective mode emerges. It is the existence of this mode 
that is the finger print [16, 23-26] of non- Fermi liquid be- 
haviour associated with the d-wave Pomeranchuk insta- 
bility in continuum and square lattice models. The emer- 
gence of z = 3 overdamped mode in our system is further 
confirmed by diagonalization of our bosonized Hamilto- 



Model Hamiltonian - We wish to describe a two-orbital 
interacting system. Hence, our starting Hamiltonian con- 
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tains a kinetic term of the form, 

(i) 

defined on a square lattice with degenerate d xz and dy Z 
orbitals per site. Here a is the spin index and r, are 
Pauli matrices. d' aa {k) creates an electron on the or- 
bital a with momentum k and spin a, and we define 

e+,-,a;y(fc) can be obtained 
by including various hopping parameters which vary from 
material to material. Since we are interested in an orbital 
ordering instability in the charge channel, only effective 
inter- and intra-orbital Coulomb interactions are consid- 
ered here. As a result, the minimal interacting Hamilto- 
nian Hi is 

Hi = Un ia ^n iai + ^ \ U' - — j n a n b , (2) 

ia i,b>a ^ 

where U and [/' are the intra- and inter-orbital interac- 
tions, and J is Hund's coupling. 

Previously, we used RPA to show that non-Landau 
damping exists in a 5-band model [16]. To set the stage 
for the bosonization calculation, we discuss briefly the 
results of an RPA analysis on the two-band model con- 
sidered here. Wc find that the self-energy of the quasi- 
particle on the Fermi surface shows a non-Fcrmi liquid 
behavior (i.e. Yl(kp,uj) ~ ui x with A < 1) in the critical 
region near the orbital ordering quantum critical point 
(OOQCP). The consistency of this result with our previ- 
ous 5-band model implies that it is the fluctuations asso- 
ciated with the d xz and d yz orbitals that leads to the non- 
Fermi liquid behaviour. Moreover, the simplicity of the 
present two-orbital model allows us to do further analysis 
on the overdamped z = 3 mode using a non-perturbative 
approach, the details of which we now present. 
Multidimensional Bosonization Multi-dimensional 
bosonization is ideally suited to this two-band problem 
because the d xz and d yz bands are quasi-ld. Following 
the standard procedure[19-22, 24], we rewrite the tight- 
binding Hamiltonian in the eigen-band index in order to 
correctly identify the Fermi surfaces and the interactions 
between quasiparticles on the Fermi surfaces. Following 
the same convention used in Rcf . [27] , we introduce a uni- 
tary matrix V av g such that the creation operators in the 

band index can be expressed as 7^ - = tl>' JJ ?, where 
v denotes a or ft Fermi surface. Using the recipe out- 
lined by Haldane, we coarse-grain the Fermi surfaces into 
N equally sized patches of width A and thickness A, as 
shown in Fig. 1. We enforce the limit of A <C A <C &f so 
that the deviation from the multidimensional bosoniza- 
tion due to the processes of momentum-transfer between 
patches and the effect of curvature within each patch can 
be significantly reduced[19]. 

In the limit of low energy and long wavelength, the en- 
ergy dispersion can be linearized near the Fermi surface, 
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FIG. 1. Illustration of the Fermi surface patches used in the mul- 
tidimensional bosonization. The tigh-binding hopping parameters 
arc ii = —1, t2 = 0.5, (3 = —0.6, t4 = —0.5, and the chemical 
potential is fi = 0.5 

effectively reducing the kinetic term to H t = v% ■ 

(k - k F )^ t1 va %- It has bcen shown[19 22, 24] that 
this Hamiltonian can be entirely described by the den- 
sity fluctuation operator defined as 

Sn s ,uq = ^Zil\, va lk+q^a ~ s 1fl n %, va )> ( 3 ) 

where the summation over momentum is restricted to be 
within patch S. Making use of the special commutation 
relation between these density fluctuation operators [28], 
we rewrite the kinetic term as 

Ht= H ^7^. Sn s^~qSn s ^ ? , (4) 

S.u,q ^ ' 

where represents summation over patches in the limit 
N — > 00 and A — > 0, which can be changed into line 
integrals along the Fermi surfaces. 

Similarly the interaction Hamiltonian in (2) can be ex- 
pressed in terms of the density fluctuation operators as 
well. After a long but straightforward calculation, we 
arrive at the normal-ordered interaction Hamiltonian, 

:Hi:=— V U*^ JJ^8nsuqU; ay m h '5n Tv _* 

a,q,ST/j,v 

(2CT_~J) y, bu y 

N k S+q k s bftg k T -q k T bU > q - 

q~STnv,a=jtb 

(5) 

where JJ^P = U~ a — c? and = —US a = sr and the 

k k k k k k 

spin index is dropped hereafter. Following Ref. [24], we 
write the effective action for the present bosonic theory, 
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S = St + Si, where 

S,va.q 

Xs,v(<f> w )= N A Q )vs ■ q/(u -v s -q) 

= N V (0) [P u VS i . + invsM- - v s .q) )(6) 



and 



/ 2lT N ^ ks+q ks k T -q k T ' 

qST^v a 

(W-J) y ^ uTu * b „ Jjbu^ S 

TV k s +q k s k T -q k T " 1 v ' 1 

(7) 

We have introduced a small imaginary part to the denom- 
inator of Xs to separate it into a real and an imaginary 
part, which will be helpful for later analysis. One can 
easily check that the interaction between quasiparticles 
with kg and kg is different from that between ks and Rkg 
(Rks denotes the new momentum obtained from rotat- 
ing ks by 7r/2). This directly means that the interactions 
contain both I = and I = 2 channels which can be de- 
coupled by introducing the auxiliary fields corresponding 
to I via Hubbard-Stratonovich transformations as 



MO) = Z) + S ks„+q S ksJ Sn ^ 



u,S 



Sv+q^kSv S ks v +q S ks, r n Sv,q: 



(8) 



where the + and — signs in A 2 correspond to the a and 
(3 band respectively and subscripts are added to patch 
labels to avoid ambiguity. Integrating out the density 
fluctuation field 5n leads to an effective action purely in 
terms of the auxiliary fields: 



M 20 M 22 l \A 2 (q) 



where 



M 22 (q,(jj) = — B + > — —^-(cr , -cr — s/ . -sr 

Z_v TV k S v +q k Sv k Sv +q k Sv 



X ^ C ks v —q"ks v S ks v —q S ks, J ^'' 



(9) 



and B = U/2 — U' + J/2. It is important to recognize 
that the A 2 (0) field is associated with the orbital ordering 
parameter which breaks the C4 symmetry. To see this, 
one can exploit the unitary matrix U c to transform 



A 2 (0) back to the orbital basis, and the resulting quantity 
will give the difference between occupation the number 
of the yz orbital and xz orbitals. As a result, we will 
focus on the region near OOQCP, that is, M 22 (0) w 0, 
and the collective modes, if any, can be determined by 
the condition M§ 2 - M 00 M 22 = 0. 

To evaluate the OOQCP condition, M 22 (0) = 0, we 
take the limit co/q — » and then q — > as advocated 
previously [24]. As a result, Rexs(O) = -N(0), and the 
condition for the OOQCP is 

^-1 + (U - 2U' + J)i ^2 N -(0)I^j > 0, 



4. r 



(10) 



Indeed, our condition for the OOQCP given in Eq. (10) 
is a generalization of the condition for the d-wave Pomer- 
anchuk instability, f 2 N(0) < —1, in the continuous 
model[24, 29]. 

Now we turn to the collective modes in the critical 
region near the OOQCP. The low energy and long wave- 
length limit corresponds to q — > and ui/{qvs) 0, 
where vs is the average Fermi velocity. A small q expan- 
sion on M 22 gives 



M 2 2(q,<- 



-B- 



■E 

v,S„ 



B 2 X°sM") 



N 



4J + 0(q*) 
(11) 



The q 2 term can be separated into real and imaginary 
part by the same trick. However, the imaginary part is 
higher order and can be neglected. Performing a similar 
analysis on Mqq and M02, one can find that in the small 
q and u>/(qvs) limit, 



M 22 
M 00 



M. 



M. 



(o) 
22 

(0) 
00 



.LO ■ 



i-M, 



(0) 



, 22 +M 2 ( 2 V + 

Q 



'00 



M ( V 



M 02 =Mg ) q 



(12) 



Consequently, we find that near the OOQCP (M$ « 0), 
the solution to 



M^ 2 )q 2 = (13) 



defines the collective z = 3 overdamped collective mode. 
This mode and has a strong dependences on the Fermi 
surface topology and momentum q. In the low-energy 
limit, oj/(qvs) 0, the condition uj / q = v$ ■ q basically 
requires that vs is perpendicular to q. Therefore, the 
Fermi surface should be smooth enough such that for an 
arbitrary direction of q, there exists at least one perpen- 
dicular vs- This is not always the case, for example when 
the Fermi surface is a perfect square. For realistic models, 
Af 22 ^ is always finite except when q x = ±q y . This obtains 
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FIG. 2. Eigenvalues of the imaginary mode of the bosonized 
Hamiltonian for BN(0) = —4.344 and q y = 0. The black dashed 
curve is a the fitting with the functional form, aq x . The case 
of q x < 0.1 is beyond the numerical accuracy with the choice of 
2000 patches. More patches are required to access to the region of 
smaller q x . 

because c- — si vanishes when fcc T = ±fc«,,. There- 

k s _ k s °' °' v 

fore, there will be no overdamped modes along the Bril- 
louin zone diagonal, which matches precisely with previ- 
ous studies of the Pomeranchuk instability on a square 
lattice [25, 26]. 

It is worth making a comparison between our re- 
sult and the previous study on a continuous model by 
Lawler, et. al. [24]. They demonstrated that the z = 3 
overdamped collective mode emerges close to the criti- 
cal point in the continuum model when an interaction 
is present in the I = 2 channel, which is similar to 
the case of the itinerant ferromagnetic quantum criti- 
cal point[30, 31]. It is remarkable to see that such an 
overdamped z = 3 collective mode exists in our lat- 
tice model as well, which strongly suggests that the or- 
bital order in a lattice model is essentially equivalent to 
the nematic order in a continuous model. Furthermore, 
the existence of this overdamped z = 3 collective mode 
from the non-perturbative multidimensional bosoniza- 
tion technique builds a solid foundation for non-Fermi 
liquid behavior since the single-particle Green function 
is changed fundamentally and obtains a non-perturbative 
form in the presence of this mode, as shown by Lawler 
et. al\2A}. 

Numerical Result - The collective modes can also be ob- 
tained by performing a generalized Bogoliubov transfor- 
mation on the bosonized Hamiltonian which allows us to 
make a direct comparison with the semi-analytical result 
obtained above. For demonstration purposes, we choose 
a set of model parameters given in Fig. 1 which have 
hole pockets a\ and ct2 only. We have just considered the 
fluctuations on the ot\ hole Fermi pocket which is suffi- 
cient to capture the emergence of the z — 3 overdamped 
collective mode. Adding up Eqs. 4 and 5, we obtain the 
resulting Hamiltonian, 



H = 



1 



N(0) 



{S s ,t + N(0)U s ,t($) Snr^nsj, (14) 

ST,q 



The density fluctuation operator Sn can be 
rewritten in terms of bosonic creation and annihilation 
operators [20, 24] 



where Us,t{<1) = B/N(c^ 



S+q°ks 



S k S +q 



l S k s )( C k s ~q C ks 



5n S -q = \/<f • v s a s ,q^[q- Vs] + y/-q- v s a\ ^[-q- v s ] 
&n s ,q = \]-q- v s a s ^\-q- v s ) + yjq- v s a\^[q- v s ] 

It can be checked that a and a' must satisfy the stan- 
dard commutation relation for bosons in order to satisfy 
the unusual commutation relation between <5?i[20, 24]. 
The Hamiltonian can now be rewritten in terms of these 
bosonic operators and diagonalizcd with a generalized 
Bogoliubov transformation. 

The diagonalization of the bosonic Hamiltonian is 
done with 2000 Fermi surface patches and the interac- 
tion parameters are set to the values for the OOQCP, 
(U - 2U' + J)N(0) = -4.334. For each momentum 
g, we diagonalize a bosonic Hamiltonian with a size of 
4000 x 4000. The energy of the overdamped collective 
mode can be identified uniquely as the only purely imag- 
inary eigenvalue of the bosonized Hamiltonian for each 
q. Fig. 2 plots the magnitude of this purely imaginary 
eigenvalue Xh as a function of q x for q y = 0, which can be 
fitted perfectly with a function of the form aq x (dashed 
curve). This proves that this branch of the overdamped 
collective modes indeed has z = 3. We have also checked 
another choice of model parameters given by Qi et. al. [27] 
as a minimal model for iron-based superconductors. In 
this case, the electron pockets have a much larger den- 
sity of states than the hole pockets, and we find that the 
OOQCP is given by U/At ps 1.7, which is in a reason- 
able range to be experimentally relevant. We still find 
the same z = 3 overdamped collective mode from the 
technique presented above, which supports our overall 
conclusion that the overdamped critcal mode with 2 = 3 
is a general feature in a two-orbital model close to the 
OOQCP. 

Conclusion - Using the non-perturbative multidimen- 
sional bosonization, we have demonstrated the emergence 
of a z = 3 overdamped collective mode from a gen- 
eral two-orbital model in the vicinity of the orbital or- 
dering quantum critical point. Since it has been well- 
established that the very existence of a z = 3 over- 
damped modes [23-26] completely washes out the stan- 
dard Fermi liquid description, a non-Fermi liquid be- 
havior should generally occur in a two-orbital model 
or in a multiorbital model with active degenerate d xz 
and d yz orbitals. Our bosonic theory provides a solid 
and non-perturbative foundation for the interpretation 
of the anomalous zero-bias enhancement observed in re- 
cent point-contact spectroscopy experiments on a variety 
of the iron-based superconductors[3, 32] as non- Fermi liq- 
uid behavior induced by orbital fluctuations[16]. 
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